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results indicate that improvements in the coloring algorithm 
may yield acceptable schedules. 
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1. Introduction. 

The procedure of scheduling classes varies with each 
academic institution. In most instances schools announce 
beforehand a list of courses and times they are to be 
offered. In this case students schedule themselves to 
courses which do not conflict. A different scheduling 
problem - occurs in the case where students are required to 
take -predetermined courses. Here the problem is to 
schedule the courses so that students have no conflicts. 

The problem of scheduling classes at the Naval Postgraduate 
School, Monterey, California, falls in this category. It 
is this problem which is investigated here. 

The digital computer has been used as an aid in 
scheduling before. Appleby, Blake, and Newman id attempt 
to produce school time-tables using a general-purpose 
digital computer. Their suggested methods of approach 
include 

(1) a random trial solution updated to eliminate 
conflicts, 

(2) a trial of all combinations possible until one 
with no conflicts is found, 

(3) a random buildup of entries into a schedule, 
adding entries only when conflicts do not occur, modifying 
previous entries to make new entries possible, and 

(4) a heuristic approach. 

The importance of knowing whether or not a solution does 
exist is stressed. The last of three successively more 
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complex programs is estimated to speed up scheduling by a 
factor of 200 to 400 over the hand operation. 

Csima and Gottlieb C6J describe a time-table construc- 
tion technique which uses matrices of zeroes and ones to 
produce results in the limited situations attempted. In 
real situations where preassignments or special constraints 
appear, further progress must be made. It is also hoped 
that with refinements in the program the computer execution 
time can be reduced. 

A method of determining examination time-tables has 
been described by Cole C5J. Citing Sherman's C9J pessimism 
in finding a schedule for a large school, he emphasizes 
producing an acceptable solution for as many courses as 
possible with reasonable computer time. Dependent on the 
imposed constraints, the execution time on the small-store 
computer was estimated to be several hours if all of 3^-0 
courses were considered. 

Another method using GASP, Generalized Academic 
Simulation Program, a program written at M.I.T. C7], 
approaches the problem by assigning students to a previously 
determined course schedule. 

This paper will treat the scheduling problem by 
assigning days and hours to courses, so that the 
constraints imposed by the students as well as instructors 
are satisfied. 
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At the Naval Postgraduate School there are aspects to 
the problem which make it different from most institutions. 
These aspects are treated as restrictions to the basic 
scheduling problem and will have an effect on its solution. 
For instance, the government allots a given length of time 
to students to follow a prescribed curriculum of study. 
Since certain courses are taught in sequence, students 
are required to take certain courses at certain times. 

The students' schedule thus loses some flexibility. 

Another aspect to be considered is the fact that all naval 
aviators must be scheduled either a complete morning or 
afternoon for weekly aircraft flights. Other restrictions 
include days off for instructor research and allowances 
made for faculty and staff education at the school itself. 

A graph theoretic approach, which makes use of the 
computer for the necessary calculations, is used to attempt 
the solution to the scheduling problem. A graph theoretic 
model of the scheduling problem is constructed. The 
problem of scheduling courses without conflicts is shown to 
be equivalent to a graph coloring problem. A digital 
computer program is produced to solve the coloring problem. 

Section 2 of this paper defines the problem investi- 
gated. Preliminary definitions and terminology are first 
presented in graph theoretic terms. The graph theoretic 
model is then constructed with these tools. In section 
3 an algorithm for graph coloring is analyzed. This 
algorithm then is used to construct a coloring program. 



9 



The results of this program are displayed in Table 1. 
Observations and conclusions based on these results are 
found in Section 4. The Appendix contains both general 
and detailed flowcharts, which explain the coloring 
program used. 
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2. Definition of the Problem. 

A. Preliminary Definitions. 

We wish to consider the basic problem of scheduling 
from a graph theoretic standpoint. But first we need the 
following basic definitions taken from Berge C 33, Chapters 
1 and 2 . 

Definition 1. Given the sets X and Y. A function iTJ 
mapping X into Y, is multi-valued , if to each xeX there 
is a well defined subset iPx^Y. Note: IP x may be the 
null set. 

Definition 2_. |Px is said to be the set of followers of 

x. 

Definition 3. A pair denoted by G=(X,)P) is said to be a 
graph whenever 

(l) X is a set, 

and (2) ^ is a multi-valued function mapping X into X. 

It is customary to represent the elements of the set X as 
points in a plane. 

Definition 4. If G=(X, |P) is a graph and x6 X, then x is 
said to be a point , vertex , or node of the graph G. 
Definition 5_. If x€ X and y€ X such that yg r* , the 
ordered pair (x,y) is said to be an arc of the graph G. 
Arcs of a graph are usually represented by an arrow from 
the first to the second element of the ordered pair. 

In case the order in which the vertices appear in the 
ordered pair (x,y) has no significance, we denote the pair 
by £x, yj=£y,x]. Such an unordered pair is called an edge 
of the graph G. 
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An edge of the graph is usually represented by a 



curve joining two vertices. 

Definition 6. If G=(X,V~’) is a graph and U= {[x,y] :y6 f 
for all xj the set of edges of G, then an alternate 
notation for the graph is the pair (X,U). Since this 
graph contains only undirected edges, it is sometimes 
referred to as an undirected graph. 

Associated with every graph is a matrix. We shall 
define it as follows: 

Definition 7» If G=(X,U) is a graph, then A=(a^^) is 
called the matrix associated with G if a^ ^ equals the 
number of edges from node x^ to node x^. 

Note: A is a square symmetric matrix. 

Definition 8. Given the graph (X,U), where x€ X and 
y6 X, x/y, x and y are said to be adjacent if [x,yj = 

[y,x] e u. 



X 



Definition 9. If G=(X,U) is a graph and if x6X, the 
number of edges having an endpoint at x is said to be 
the degree of x. 
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B. Definition of the Scheduling Problem. 



Definition 10 . A section is a group of students all of 
whom take the same courses. 

Sections are used to simplify the scheduling problem. 
Determining a single section schedule is easier than 
finding a separate schedule for each individual. 

Denote sections as s. , i=l,...,n,, courses as c 

J 

j=l,...,n 2 * and professors as p^, k^l^...^^. We shall 
consider the set X=|s^J \J j as a set of nodes for 

a graph. We shall also consider the sets of all pairs 
U^= , c jJ : section s^ must take course c^J and 

p^: instructor p^. is assigned to teach course c^.J . 

Define U=U^UUg. We see that U with the set X defines 
the graph (X,U). 

Definition 11 . The grapn G=(X,U) described above is 
called the data graph . 

We shall consider an example with the following 
conditions : 



(1) Section s^ is required to take courses c^Cg, 
and c^. 

(2) Section Sg is required to take courses c^ and 

c 2| • 

( 3 ) Instructor p^ is assigned to teach courses 
c-^ and c^. 

(4) Instructor p 2 is assigned to teach courses 
c 2 a nd c j j . 
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Then sis deTineci o/bovs^ X — 

u= ^[?i> ^1 3 c 2 ^ } C s i -* c 3^ -» E s 2 1 c i±3’ C s 2 ’ c 3T rPr c iT 

[Pl' c 3 ]> I?2 ■* °2-l ’ [?2 ,C h The data graph G=(X,U) is 
illustrated as follows in Figure 1. 

P, 

Pz 

Gr 

Fig. 1. Example Data Graph 

Definition 12 . Given the data graph G=(X,U) described 
above, a conflict is said to occur between classes c^ 
and c y i^j, 

if (l) [s k , c i] €U and rv«; € U for some k, 
or if (2) C c i J Pi-J^ u and [c^pJ^U for some 1. 

This is the same as saying these classes c^ and c^ 
cannot be scheduled simultaneously because: 

(1) the same section s^ must take both classes and 
could not be present at two places simultaneously, or 

(2) the same Instructor p^ must teach both classes 
and could not do so simultaneously . 




We consider now the concept of a conflict graph 
deduced from the data graph. 

Definition 13 . Given the data graph (X,U), a new graph 
(C,U') is called a conflict graph if C=-£c^:c^ is a node 
corresponding to class c^ for all i J’ and U , = : 

classes c^ and c^ have a conflict J. Note: c^ cannot have 
a conflict with itself; therefore, for all i, E°i 3 c ±l^ U ' 
The associated matrix with the conflict graph we 
shall call the conflict matrix. The conflict matrix 
M=(m 1 j), i, j=l, . . . ,rV), is such that 

m^j= fo, if c^ and c^ have no conflict. 

| 1, if c^ and c^ have a conflict. 

If we consider the same example as before, we can 
display the conflict graph in Figure 2 below. 




Cr 



Fig. 2. Example Conflict Graph 

The conflict matrix in this case is M= lot /o 

1 / 0 / / 

\'J 0/ 

\o ( / o t 
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Coloring a graph is a process of assigning colors to 



each node. 

Definition 14 . Given an integer p, a graph G is said to 
he p -chromatic , if the nodes can be colored with p 
distinct colors so that no two adjacent nodes are assigned 
the same color. 

Definition 15 . The smallest number p for which a graph 
is p-chromatic is called the chromatic numb er of the 
graph G and is written k(G). 

Given a group of classes with arbitrary conflicts 
among them, we solve the problem of scheduling by 
assigning time periods to these classes. Where a 
conflict arises between two classes we require that they 
be assigned different time periods. On the other hand, 
if we consider these classes as nodes of a graph and 
their conflicts as edges, the problem of coloring the 
nodes of this graph is seen to be equivalent to the 
above assignment of time periods. 

The problem posed here is not necessarily one of 
determining the chromatic number of a conflict graph. 

Our problem is to find a suitable coloring of the conflict 
graph, which is consistent with imposed restrictions 
which are now discussed. 

The Naval Postgraduate School uses the quarter system 
for teaching courses. Each course lasts twelve weeks. 
During the twelfth week of each quarter final exams are 
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held. The first eleven weeks have the same schedule. 
Normally each course will meet several hours each week 
depending upon the nature of the course. 

Earlier we assigned a node c^ to each course. We 
shall now introduce nodes with two subscripts to 
differentiate between the hours of a given course. 

Consider a three hour course c^„ Then the hours or 
nodes associated with these hours are denoted c. ., 

1 ^ X 

c. 0 , c. oi or simply c, ., j=l,2,3. If a course is 

a one hour course, we shall let j=0. 

In some cases, if there is a large number of 

students taking the same course in a given quarter, 

it may be necessary to divide this group into two or 

more subgroups for teaching purposes. These subgroups 

are assigned to take different sessions of the same 

course. These different sessions we shall call segments 

of that course. We shall introduce nodes with three 

subscripts to differentiate between segments of a course. 

In the example above, if course c^ has two segments, then 

the hourly nodes associated with the course are c. . , , 

1 3 J )£■ 

i denoting the course, j being the hour, and k identifying 
the segment. In this case they are c. . , , j=l,2,3, 
k=l,2. If a course only has one segment we shall let 
k=0 . 

Now we let X= [ s i]^{ c j,k,lju( p m] for a11 k > 

1, and m. We also consider U^=^|j3^, c ^ ^ J : section s^ 

must attend hour c. , .f and U 0 =ifc. , . ,p 1: instructor 

J,k,lj 2 | j , k , 1 m J 
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is assigned to teach hour c^. ^ xj* ^ U=U^ ,, then 
we have G=(X,U) as our data graph. We assume that as the 
need for segments arises, the appropriate section to 
segment and segment to professor assignments are made 
by the academic departments. We may thereafter regard 
these segments as different courses. 

Given our new data graph, we define conflicts in 
exactly the same manner as in Definition 12, except 
that we now allow course nodes to have three subscripts. 

We now consider the previous example with the 
modifications below. 

(1) Let s^ take courses c^, Cg, and the first 
segment of course c^. 

(2) Let s 2 take course c^ and the second segment 
of course c^. 

(3) Courses c^, Cg , and c^ will all be three hour 
recitation courses. 

(4) Course c^ will consist of three recitation 
hours and two laboratory hours to be scheduled together. 

(5) Instructor p^ will teach both segments of c . 



Department meetings, instructor research time, 
aviator flight time, and other such periods imposing 
constraints are displayed on the conflict graphs. In 




l,i,0' c 2,i,0' c 3,l,j ,c 4,k,0' 



Figure 3 illustrates G, the data 



graph. 
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Pig. 3. Data Graoh 
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the case of a one hour lecture for all students, we add 



to X a lecture node c, n n for some 1, and to U we add 



edges £ c 1 Q q> s iJ for all 1 • In the case of a one hour 



meeting for all professors in a given department, we 
prM -t-n y o nnHo r> fnr> some k , and to U we add the 



>ly other constraints are 

represented . 

Our data graph is G=(X,U). X is the set of all 
nodes denoting sections, professors, courses and their 
segments, flights, instructor research, and special hours 
such as department meetings or special lectures. The 
set U consists of all applicable edges. 

Using the above data we must introduce a slightly 
different concept. 

Definition 1$ . Given the data graph G=(X,U), a new graph 

G'=(C,U') is called an hourly conflict graph for G, if 

C=/c ± ^j c i j k is a node corresponding to course c ±J 

hour j, and segment k } and U' = /(c, , , ,c. ~J: hours 

J ^ 1 j j jK- ± jiiij n 

c. . , and c, m have a conflict} 
i,j,k l,m,n J 

Based on our previous example the hourly conflict 
graph is illustrated in Figure 4. 




instructor p^ belongs to 
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Pig. 4. Hourly Conflict Graph 



Given a course that meets several times a week, we 
consider now conflicts between the hours of the course 
described above. These hours may be divided into 
recitations or laboratory work- For a given course, 
using the second subscript, we shall number the recitation 
hours first. For instance, a non-segmented course c^, 
a five hour course consisting of three recitations, will 
have its hours denoted c. . _, c. Q _, c. - n , c. h _, 
and c. n , the first three of which correspond to 
recitations. Note: [c. . A ,c. , A l€U' for all jj^k. 

L 1 j JjU 1 , K , U J ’ 
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We shall restrict our problem somewhat by allowing 
only one recitation per course to be held daily. This 
brings about the following definition: 

Definition 17 . If C is the set of all course hour nodes, 
then there exists a daily conflict between c i ^ and 



and (4) j and m both correspond to recitation hours. 

We now augment our information regarding conflicts 
by the following: 

Definition 18 . If G=(X,U) is the data graph, then 
G' '=(C,U' 1 ) is called the daily conflict graph for G, if 
C is the set of all course hour nodes and U''=<fc. , , , 



If we consider our previous example again, we see 
that Figure 5 below illustrates G 1 ', the daily conflict 
graph. 



c. provided: 

l,m,n y 



( 1 ) i=l, 

(2) k=n, 

(3) #m. 




have a 



daily conflict 




Fig. 5 . Daily Conflict Graph 
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In the case of laboratory hours of a course or flight 
hours for aviators , these hours must be scheduled together 
during a given day. To consider this we introduce the 
following : 

Definition 19 . If G=(X,U) is the data graph, then S= 

(C,U' 1 1 ) is called the sequence graph for G, if C is the 

set of all course hour nodes and U ,,, = -{[c. . , ,c, m " 7 : 

( L i Jjk l,m,n J 

c. . , must be scheduled together with c, v. 

The sequence graph S for our previous example is 
shown in Figure 6 below. 

( “W Q-i 2 , 0 ^4,5,0 

# # • # 

4» 4° 4v 9v c v< j 4° 4° 

Fig. 6. Sequence Graph 
We have a data graph which summarizes our input 
information containing all section, course, and instructor 
assignments. The hourly conflict graph G' displays course 
hours which must be scheduled on different hours. The 
daily conflict graph G' 1 displays course hours which must 
be scheduled on different days,. Our sequence graph S 
indicates which course hours must be scheduled together 
in block form. The problem is to schedule these course 
hours in a five day week. The school allows nine class 
hours per day. 
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We will now consider the concept of a color vector 
with two components, the first of which will be a color 
corresponding to the day. The second will correspond to 
the hour within the day. Our problem is to color the 
course nodes of the hourly conflict graph G' so that no 
two adjacent course nodes have the same color vector. 

The coloring of the nodes of G' may be thought of as a 
simultaneous coloring of all three graphs G 1 , G 11 , and 
S subject to the following restrictions: 

(l) no two adjacent nodes of G 1 are assigned the 
same color vector, 

and (2) no two adjacent nodes of G 1 1 are assigned the 
same day color, 

and ( 3 ) adjacent nodes of S are assigned the same day 
colors and sequential hour colors. 

Note: It is sufficient for two color vectors to be 

different when either or both of the corresponding 
components are different. In view of this it is interesting 
to note that satisfying (2) above automatically satisfies 
(l) for the same two nodes. 

In this coloring process we are restricted to five 
day colors since the schedule desired is a weekly 
schedule. In addition, the number of hour colors within 
a day must be restricted to nine. Given a satisfactory 
coloring of the conflict graph, we assign days to the 
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day colors and hours to the hour colors. The resulting 
schedule meets all the constraints Imposed provided 
we have not exceeded the maximum number of colors 
allowable In each case. 
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3. The Coloring Process. 

A, The Algorithm, 

Although we cannot be certain of obtaining the 
chromatic number for a given arbitrary graph, coloration 
solutions may be attained by trial and error methods or by 
linear programming techniques f 33. An algorithm for 
coloring a graph, introduced by Welsh and Powell 03 , has 
been used to attempt to solve the vector color problem of 
the Naval Postgraduate School course graph. This 
algorithm is described as follows: 

Let G=(X,U) be an arbitrary graph where X is the set 
of nodes and U the set of edges. Denote the vertices X= 
■^Aj > i=l,...,nj. Let the degree of the vertex A^ be 
denoted by d^, where i=l,.,.,n. Assume that the vertices 
have been numbered so that their degrees form a non- 
increasing sequence d^ £ d g ^ , , , ^d (j , We now define T.CX, 
where i=l,,,,,t, T^ is defined as follows: 



( 1 ) 




*■ k**i 

for some m, then A.€T 

J 



1^, where i^l and i < i ? < . . ,<i 
1 , if and only if, 



(a) j>i , 



m 



and (b) [A ,A. } d U for k-l,.,.,m. 

J 9 

Note: T^ 0, 

Similarly we define T g , where T 2 C(X-T^). is 



constructed as follows: 



(l) If i is the least integer for which A.^ T^ , 
then A^€ T^, 
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(2) If J^e Tg, where J 1 < J 2 < , . ,<J p , then 

A^€Tg, if and only if, 

(a) /> j. D , 



and 



0>) f A * ' a j, 3 i U where k=l,...,p. 



If T p =0, we are finished. If not, we construct T, C (X— 

l ^ 3 

C/T, ) in the same way as we constructed T p . 

i 1 ^ 

Since X is finite, there will exist a least integer 

r-/ 

r such that T r =0. We note also that tjT^X and that if 
i^j, then T^f)Tj=0. We can now assign a different color 
to each If . If A^,Aj€T k , i^j, for some k, then Qa^,A 
<j^U. Thus no adjacent nodes will have the same color. 

The algorithm furnishes an upper bound for the chromatic 
number k(G). Thus k(G)^r-l. 

Welsh and Powell further state that 
k(G)-^ r-l^max min(i,d.+l) 

i 1 

Proof of this upper bound comes from analysis of the 
algorithm. By construction, each node belonging to 
(X-T^) is adjacent to a node of T^ # By construction 

Ml 

we see that each node of (X-^T.) is adjacent to at least 



one 



d k ^ m. 



f 7 Ml 1 *** 

node of each of | T^.j , Hence, A k 6(X-C/T^) implies 



d;+l 



Consider (X-^/T.). A. cannot belong to this set 
l J 1 



because d. >d.+l. Hence, 
L d;tl 

(1) J.6UT,. 



Since each T. contains at least one element, 

J » 

(2) A6L/T 

I J 
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Together (l) and (2) above imply that 

W) 

(3) A. e t/T., where 0?(i )=mln (i,d,+l ) . 

1 • j c/cj - c m 

Let o (G)=roax (0(i) . Then A. € L/T l^T for all i . 

i J 1 i J I J 

Setting q= <f (g)+i we see that G^=(2f. Thus for the least 
integer r where T=(2f, this algorithm ensures r^q. Hence^ 

k(G)^r-l ^q-l= ^"(G)=m^.x min(i,d 1 +l). The desired 
inequality is proved. 

The above 6 (G) provides us with an upper bound for 
the chromatic number of a graph. Berge's [2] statement of 
Brooks ' 01 theorem on coloring the nodes of a network 
gives us insight about a lower bound. This statement, 
simplified to apply to our problem, is as follows: 

Theorem (Brooks ) . A graph G of maximum degree q where 
q>2 can be colored with q colors except in the case where 
a component of G consists of q-fl vertices, each of which 
is joined to the other q vertices. 

The above mentioned component is referred to as a 
complete subgraph of G. In other words, if G contains a 
complete subgraph on m vertices, then k(G)^m. Thus the 
largest complete subgraph of G sets a lower bound on k(G). 

To illustrate the use of this algorithm in coloring, 
consider the given arbitrary graph in Figure 7 below. 




Fig. 7. An Arbitrary Graph 
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Numbering of the nodes must be done so that their 



degrees form a nonincreasing sequence. One such 
renumbering is listed below. 



Node 



Degree Renumbered 



2 



6 



4 



1 



5 



3 



4 



4 



3 



3 



3 



3 




Figure 8 below illustrates the arbitrary graph with 



Az 

Fig. 8. Relabeled Arbitrary Graph 
Using the algorithm to form the disjoint sets lb, 
i=l„ > 4 ) we see the following: 



To obtain a satisfactory coloring, we assign the nodes of 
T^, Tg, and colors 1, 2, and 3 respectively. 

We note from the above that this arbitrary graph is 
3 colorable by the algorithm. Of interest is a comparison 
of this number with the lower and upper bounds of the 
chromatic number for this graph. By observation, the 



nodes relabeled. 




T 



1 



|Ai,A 2 |, T 2 |a 3 ,A 5 |, 



, and T^=0. 
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largest complete subgraph has 3 nodes; therefore, the lower 
bound is 3. Since the number of colors in our solution 
equals the lower bound, we deduce that the chromatic 
number is in fact 3. The upper bound of k(G), cf (G)=max 
min (ijd^+l )=max|l,2 , 3* 4, 4 *4 J - is 4 for this example. 

While in general we will not know what k(G) is for an 
arbitrary graph, this algorithm does ensure a p-coloring 
where p<^<^(G)=max min (i,d^+l). 



V 
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B. Coloring the Graph by Color Vectors. 

In coloring our conflict and sequence graphs we are 
restricted to five day colors and nine hour colors per 
day. Flowchart A in the Appendix illustrates the coloring 
procedure and application of the Welsh and Powell algorithm 
to this problem by means of a digital computer program 
using Fortran IV language. This procedure is described 
briefly in chronological order as follows: 

(1) Section and instructor course assignments are 
read in. 

(2) A course directory is established. 

( 3 ) An indexed list of conflicts is determined based 
on section and instructor conflicts, 

(4) Course hour data is read in, 

( 5 ) Course hour nodes are given day colors using an 
arbitrary selection process. 

(6) Course hour nodes for each day are given hour 
colors using the algorithm. 

Data for the program has three forms. First are the 
sections cards containing the section, the number of 
students, the number of aviators for that section, and 
each course taken by that section. To differentiate 
between courses, segmentation and lab designations were 
used. To indicate instructor assignments, data cards were 
used with the same format as on section course assignment 
cards. This procedure simplified the computation of 
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conflicts between course hour nodes. Flowchart B in the 
Appendix illustrates the read-in and print-out of conflict 
data. 

Flowchart C in the Appendix Illustrates the 
computation of a course directory. Each course is 
associated with an integer, and the integer listing of 
these courses is used for all further computations. This 
course listing is printed out for reference. Included in 
this computation is a random number generator used for 
random reordering of this course directory. Flowchart 
I in the Appendix illustrates RAN the generator function 
used . 

The list of conflicts is computed next. This is 
shown in Flowchart D in the Appendix, In the example 
described in Section 2, the directory assumes the 
following form. 

CC(I) I 

C0URS1 1 

C0URS2 2 

C0URS31 3 

C0URS32 4 

C0URS4 5 

C0URS4L 6 

C0URS3I and COURS32 are segments 1 and 2 respectively of 
C OUR S3 and C0URS4L is the lab portion of C0URS4. The 
conflict list for the above example is shown below. 
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J 



1234567 

INDEX (J ) 1 4 8 11 15 18 21 

I 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 ... 

LIST (I) 234135612413562... 

The above information is interpreted as follows: 

C0URS1 whose integer representation is 1 has conflicts 
with C0URS2, C0URS31, and C0URS32, whose integer 
representations are 2, 3 > and 4 respectively. The 
conflicts are found by entering the list at index (l) and 
reading all entries up to index(2), where the conflicting 
courses of C0URS2 begin. This notation has the advantage 
of requiring less memory space and less print-out space 
than the comparable matrix representation discussed 
earlier. The nonzero elements of the I-th row or 
column of this matrix are listed starting at index(l). 

The list which describes the conflicts between the 
nodes is ordered so that the nonzero elements corresponding 
to the I-th row or column of the matrix form an increasing 
sequence. Flowchart E in the Appendix describes this 
procedure and the print -out of its results. Flowchart 
J in the Appendix displays the subroutine RELIST used in 
this procedure. 

Incidental to the main objective of the program, an 
upper bound for the chromatic number of the course graph 
is computed. During the development of the program it 
served to give an indication of the degree of difficulty 
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to be encountered later in the coloring phase. Flowchart 
F in the Appendix shows the computation of the upper bound 
for this graph using courses' as nodes. 

We next read course hour information into our program, 
in order that day colors may be assigned to the course 
hour nodes. Flowchart G in the Appendix illustrates this 
input and assignment of day colors. Flowchart K describes 
the subroutine DAYCOL called for assignment of day colors. 
This routine assigns recitation hours of the same course 
different day colors. Each lab course is recognized as 
such in order that its hours may be scheduled during the 
same day. A listing with index is compiled as shown, to 
remember these day color assignments. In our example the 
course hour information is as follows: 



I 


CC(I) 


HOURS 


1 


C0URS1 


3 


recitation 


2 


C0URS2 


3 


recitation 


3 


C0URS31 


3 


recitation 


4 


C0URS32 


3 


recitation 


5 


C0URS4 


3 


recitation 


6 


C0URS4L 


2 


lab 


Using the above 


course hour 


information DAYCOL 



assigned the day colors in our example as follows: 
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To each day we now apply the coloring algorithm. 
Flowchart H in the Appendix gives the general procedure 
for the hour coloring. Flowchart L illustrates the 
HRCOLR subroutine used in that procedure. When labs 
having more than one period are scheduled on any given 
day., the integer representation for that course will 
be repeated for that day. To enable these nodes to 
maintain their separate identities, a renumbering of 
nodes is made. Relative to this new numbering, an indexed 
list of conflicts is computed from the master list for 
each particular day. Then the degree of each node is 
determined. Knowing this, the program assigns colors to 
each node according to the coloring procedure described 
earlier. The result is the second or hour component of 
a color vector. Having day colors already determined by 
DAYCOL, the process is complete. 

Using the above program, data from the second 
quarter of the 1967-1968 academic year at the Naval Post- 
graduate School was processed. No special periods such as 
flight hours, instructor research time, or departmental 
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meetings were included in this data. Included were 424 
courses, 403 sections, and 194 instructors. The number 
of conflicts for any course varied from 0 to 52 ; however, 
over half of the courses had conflicts ranging from 5 to 
13. 

The number of hour colors required for the five days 
were 9 > 8 , 10, 10, and 10, To investigate the effect of 
the DAYCOL assignments on the total number of hour colors, 
a random number generator was used to reorder the 
directory of courses prior to the DAYCOL routine in the 
program. Flowchart J in the Appendix illustrates this 
procedure for the given input data set. Table 1 
displays the results of that investigation for 18 runs, 
including the first already mentioned. Each entry is 
the total number of hour colors for a particular run 
and day color. 

None of the runs, which averaged about 7 minutes of 
IBM 360/67 time, gives the acceptable nine or less hour 
colors for all days. Column 6 of the table indicates the 
number of acceptable days for each run. Runs 3 and 18 
were the best in this category having 3 each. Runs 1 , 

15, and 17 had 2 acceptable days each. All others had 
one or less. Column 7 of the table indicates the number 
of hour colors exceeding 9 each day, summed over the five 
days for each run. 
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4. Conclusions. 

While no acceptable coloring solution has been found 
directly from the described program,, several program runs 
are sufficiently close to indicate that modifications in 
the coloring technique may in fact give a satisfactory 
solution where an acceptable trial and error solution is 
known to exist. The known solution, hand generated from 
the same data with all constraints taken into account, 
used nine periods for each of the five days. Coloring runs 
1 and 18, Table 1, both came within 3 periods or hour 
colors of equalling the number of periods required within 
one week. Huns 3# 12, and 17 were also noticeably close 
to a solution. 

In this case where a solution is known to exist, the 
results obtained must be attributed to the coloring 
technique applied. With an acceptable solution in mind, 
we see that any coloring algorithm or technique which yields 
nine or less hour colors per day is a satisfactory pro- 
cedure. Changes in the present program in DAYCOL or in 
the algorithm itself could improve results. One such 
change in the algorithm may be the following: 

T^ can be determined by the Welsh and Powell algorithm. 
The set of nodes in T^ is then removed from the graph, 
together with the associated edges. The remaining sub- 
graph is treated as a new graph, and the Welsh and Powell 
is reapplied to find the set Tg, etc. 
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Also a change, which should be considered, is a 
coloring program where the day and hour colors, dependent 
upon one another, are assigned together. This involves 
combining the functions of DAYCOL and HRCOLR, with whatr- 
ever algorithm used, into a program operating on a total 
conflict graph rather than separate hourly and daily 
graphs. 

Another observation concerning our schedule problem 
can be made. Given arbitrary data, we may question the 
existence of any schedule solution. A solution to the 
presently unsolved problem of determining the chromatic 
number for an arbitrary graph will indicate if a schedule 
solution is in fact possible. 
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TABLE 1 



Day Color 


1 


2 


3 


4 


5 


Col 

6 


Col 

7 


Runs 1 


9 


8 


10 


10 


10 


2 


3 


2 


11 


10 


10 


11 


10 


0 


7 


3 


9 


9 


9 


11 


11 


3 


4 


4 


11 


9 


11 


10 


10 


1 


6 


5 


10 


10 


12 


11 


10 


0 


8 


6 


9 


12 


10 


12 


10 


1 


8 


7 


9 


13 


12 


11 


10 


1 


10 


8 


12 


10 


10 


10 


11 


0 


8 


9 


10 


8 


13 


11 


10 


1 


8 


10 


10 


10 


10 


12 


10 


0 


7 


11 


10 


9 


12 


10 


10 


1 


6 


12 


10 


9 


10 


11 


10 


1 


5 


13 


11 


10 


11 


12 


9 


1 


8 


14 


11 


12 


11 


9 


11 


1 


9 


15 


9 


14 


11 


9 


10 


2 


8 


16 


12 


11 


12 


9 


11 


1 


10 


17 


11 


9 


11 


10 


9 


2 


5 


18 


10 


9 


9 


11 


9 


3 


3 
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APPENDIX 



GENERAL COLORING PROCEDURE 
FLOWCHART A 




f 
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READING CONFLICT DATA 
FLOWCHART B 
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COMPUTING COUKSa LlGT 
FLOWCHART C 







XKQsR(l) 




YKQsCC(I) 




DO 


106 


1*1, IC2 




R(I)=XMAX 




CC(I)rCC(JK) 


»■ 


DO 


104 


J*1,IC3 




R(JK)2XKQ 




CC(JK)5YKQ 




DO 


102 


K=1.9 
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Write ( 6.12) ic ? 

•write ( 6,16 ) ( (KX(I, j ) , J* l,9) y 
16 FORMAT (5X,6I9)/I-1,IC3 



WRITE (6,32) (I ,CC(I) , 
Isl,IC2) 32 FORMAT 
( 5X.I5t2X t A7) 
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COMPUTING CONFLICT LISTING 





DO 24 I 1,IC2P1,6 

lie 1+5 

* WRITE (6.23) ( K , INDEX ( K ) , Ks.I , II ) 
23 FORMAT (12(4X,I4)) 




DO 25 I l.MF ,6 
* lie 1+5 

WRITE (6,23) (K,LIST(K) ,<(31,11) 
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REORDERING CONFLICT LIST 



FLOWCHART E 
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COMPUTING UPPER BOUND 



FLOWCHART F 





INDEX (IC 2*1 )sM 




DO 51 Isl,IC2 




DO kO 1*1. IC2 

40 KD ( I ) ? INDEX ( I+l )- INDEX ( I ) 




WRITE (6,12) 
KD( I) 
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READING IN COURSE HOURS AND COMPUTING DnY COLORS 



FLOWCHART G 





WRITE (6,23) h 

( K , INDEXD ( K ) ,K5l t 6)| 



DO SO Is 1 NTOTC , 6 
IP 55 1+5 
WRITE (6.23) 
(K,LISTD(K) ,KgI»IP5) 
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COMPUTING HOUR COLORS 
FLOW GRAPH H 



CALL HRCOLR (LIST, INDEX 
INDEXD ,LISTD , INDEXH , 
LISTH.NCOL) 






WRITE (6,8) 
(NCOL(I), 121,5) 

8 FORMAT (5(51.15)) 



mm (6.23) 

'( K , INDEXH ( K ) ,Ksl,6) 



DO 85 I = 1 NTOTC ,6 
IP5*I+5 
WRITE ( 6 , 23 ) 
(K t LISTH(K) ,KsI > IP 5 ) 
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FUNCTION RAN (I) 
FLOWCHART I 




51 



SUBROUTINE RELIST 



FLOWCHART J 





KQsLIST(J-l) 
LIST(J-l)sMIN 
LIST( JK )s KQ 




RETURN 

END 
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SUBROUTINE DAYCOL 



FLOW GRAPH K 





IN = 

MC( J,I) 




DO 14 Ksl.IN 
ND (LD)- ND (LD )+l 
MA(LD,ND(LD) )*I 







Hsl 

INDEXD(l) 




30 10 1*1.5 

INDEXD ( 1+1)5 INDEXD ( I )+ND ( I ) 


=1 




IN5ND(I) 






DO 10 Jal,IN 

LISTD(M)3MA(I,J 

M*M+1 




♦ 



NT0TC*M-1 

RETURN 

END 
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SUBROUTINE HRCOLR 



FLOWCHART L 



SUBROUTINE HRCOLR ( INDEX , LIST , INDEXD , LISTD , NCOL , IN*. 
DEXH.LISTH) DIMENSION INDEI, LIST, INDEXD, LISTD .INDEXH, 
LISTH , NN , INDEXC , LISTC , KD , KK , NZZ , NCOL , LISTT , NN2 1 NNZZ 








L=1 




INDEXC (K)sL 




DO 26 




DO 2S Kal.M 




IMl* INDEX (NN(K)) 

INIs INDEX( NN ( K )♦!)-! 


— ► 


Us l.M 





mm (^. 32 ) I 

I,M 32 FORMAT I *- 

(//7I,I3.5X.I3)| 
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DO 22 J 




JP5rJ4-5 





WRITE (6,23) (K, 
INDE.X.C (K) ,K- J , JP 5 ) 

23 FORMAT (12(4X f l4) 




22 



m 

L-l 





DO 20 J 




WRITE (6.23) 




DO IB 


*■ 


1,LM1,6 

J?5:J+5 




(K ,LISTC (K ) , 
K*J,JP5) 


>/20) A- 


J»1,M 





■ IMAX»KD(K) 

mSlMSJL 




14 



KQsKD(J) 

» KD ( J )s MAX 
KD(JK)=KQ 





kWKT'J) 




DO 35 




kk(j)skk(jk; 

KK ( JK )s KQ 


K2) * 


JP5?J+5 WRITE (6,23; 
(K t KK(K) t Ksj,JP5) 




35 



m 70 j«i;i!OT 

-HNZ(J)si 
DO 70 K" 1.6000 



4NZZ(J,K)J 
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LISTT(KK(J) )=NCO 




NZ«NZ+1 




NZZ(NZ)xJ 


► 


NNZ(NCO)»NNZ(NCO)«#*1 




NNZZ ( NCO , NNZ ( NCO ) )r KK ( J ) 




1M1= 1NDEA.C (KK( J ) ) 



IN1=INUEXC (KK( .1 1 

+ 1)-1 

umm 



JfES 






'J.EQ. 



./NO 


DU 6 J1=JP1,M 






HZMlxNZ-1 






JEST 

^KlalMl, 
f INL 



12- INDEXC ( KK ( J 1 ) ) 
H23lNDEXC(KK(Jl) 

* 1)-1 






IN 2 =NNZ(NC 0)-1 




NO *" 


DO SO K 2 =I 2 ,N 2 
DO go J 2 = 1 .IN 2 





m listtu^ijdi-nwi 

NNZ(NZ)- Jl 
NZxNZ+1 



NNZZ(NCO,NZZ 
(NCO) )-KK( Jl ) 

NNZ ( NCO )r. NNZ ( NCO }+l 


-0* 


NCO i 

NCCH 


*(io)-> 


NCO(I)=NCO-l 
LX/* 1 

DO 50 J-IM.IN 





LISTH(J) 

LISTT(LL) 




DO UO J- 1,600 
IND 8 XC( J)sO 


-*(ZoW 


bo 60 J- 1 , 

6000 




30 

RETURN 




LL-LL +1 




KK(.T )20 




60 LISTC(J )=0 




END 
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